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Abstract. 

In this paper a reduced set of the partial differential wave equations valid in the 
conversion layer close to O-mode cutoff surface and accounting for the magnetic field 
2D inhomogeneity with no restriction to an angle between the toroidal direction and 
the magnetic field direction is derived. An integral representation of a solution to the 
derived set of equation is given. For the particular case of small angle between O mode 
cut-off surface and X mode cut-off surface an explicit expressions for both the electric 
field components and the conversion coefficients are obtained and its properties are 
considered in details. 
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1. Introduction 

Nowadays the electron Bernstein waves (EBWs) having no density cut-offs and 
effectively damped even at high electron cyclotron harmonics are considered as the 
most promising candidate to provide an auxiliary heating and current drive in a dense 
plasma of a spherical tokamaks and stellarators |H EJ Ej • The EBW could be excited 
via direct conversion of X mode to Bernstein mode in a vicinity of the upper hybrid 
resonance (UHR) or via so-called O-X-B scheme. The efficiency of last scheme, as it was 
demonstrated theoretically OEHH] m the frame of ID slab model, is determined by the 
efficiency of O to X mode conversion, which can reach 100 percent value at the certain 
parallel refractive index being constant in slab geometry. Due to, in a real spherical 
tokamak's configurations, where the poloidal inhomogeneity of the magnetic field is 
important, both the parallel refractive index is no longer constant and the components 
of the dielectric tensor are functions of two co-ordinates, an analysis of the full-wave 
equations in the frame of 2D model is important. The first attempt to consider 2D 
model of OX conversion has been undertaken a couple of years ago [7j. The main 
conclusion provided by author of [Zj concerning the absence of the O mode reflection 
from the O-mode cut-off surface seems to be quite doubtful and lacks support from the 
last two papers devoted to this topic [HI E] • Unfortunately, the OX mode conversion 
in [B] has been considered in the frame of the oversimplified model ignoring as it 
does both the poloidal magnetic field and, as result, varying of the parallel refractive 
index on the magnetic surface. These effects were taken into account in jHj, where the 
the explicit expressions for the OX and XO conversion coefficients have been obtained 
simultaneously with jS] but for (p = arctan Bq/B z <C it (Bq and B z are the poloidal 
and toroidal components of the magnetic field and tp is an angle between the toroidal 
direction and the magnetic field). 

Unlike the assumption used in (p <C 7r , in spherical tokamak's plasma if is no 
longer small as it is demonstrated in figure 1 for a typical MAST tokamak discharge. 
Therefore, straightforward using the results of [3], without mentioning paper (HJ, to 
describe OX conversion in spherical tokamaks seems to be overhasty. Because of the 
intensive efforts are paid to provide an auxiliary heating in spherical tokamaks using 
OXB scheme, calculation of the OX conversion coefficient in the realistic 2D model is of 
great interest. In this paper a reduced set of the partial differential wave equations 
valid in a vicinity of the O-mode cutoff and accounting for the magnetic field 2D 
inhomogeneity with no restriction to an angle between the toroidal direction and the 
magnetic field are derived. A solution to the reduced set of equations has been obtained 
and it properties are considered in details. 

2. Physical model 

There are three effects that remain beyond the scope of the present paper. We neglect, 
first, the curvature of the magnetic field line at the magnetic surface because of its 
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local radius Rf is considerably greater than the beam radius p, second, the curvature 
of the magnetic flux surfaces assuming high localization of the conversion region, third, 
the magnetic field shear which is not important for OX conversion ^Oj- We restrict 
ourselves to the case of not extremely strong plasma density inhomogeneity L n when 
geometrical optics can be applied except possibly near cut-off layer or possibly mode 
conversion layer. One introduces two Cartesian co-ordinate systems (x, y, z) and (x, £, £) 
with distances scaled in the units c/uj and their origin located at the O-mode cut-off 
surface (figure 2). The co-ordinates x, y and z imitate the flux surface label, the poloidal 
and the toroidal co-ordinates, respectively, and axes and are along the magnetic 
field and perpendicular to it on the O-mode cut-off surface, respectively (figure 3). The 
transformation from y, z components to (, £ components convenient to represent the 
Maxwellian equations is 

y = cos(<p)C + sin(y3)£ 

z = - sin(v9)C + cos(y?)£. (1) 
One introduces the set of the wave equations for a monochromatic wave ~ exp (iut) as 
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: (E x ± iE c )/y/2, e± = l-v/(l±q),ri = l-v,v = uJ 2 p Ju\ q 
are electron Langmuir and cyclotron frequencies, respectively. 
Since a domain in a vicinity of O-mode cut off surface beyond the tokamak mid- 
plane is of interest, one can expand the plasma parameters, namely density n(x) and 
magnetic field modulus B(x, y), into the Taylor serious at r 

n(x) ~ n (1 + x/L ri ) 

B(x, y) ~ B (l + x/L bx + y/Lby), 

where 

l/L n ,p/L h y < 1, 

L~ l = dlnn e /dx\ ro ,L^ = d\nB /dy\ ro , 
are parameters being the first order quantity (0(1)) and 

l/L bx < 1, (5) 
L bx = dlnB/dx\ ro 

is one being the second order quantity (0(2)). In order to study the properties of the 
waves in the mode conversion region we keep in mind that the component E_ is small 
compared to two others E- ~ (E + ,E^) / L n there. Omitting in (2) the terms being 
higher order quantity than the first one obtains 
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Due to tokamak symmetry we may assume that the wave fields vary as exp (in z z), where 
n z being constant is large. We would like to construct a solution E = (E + ,E^) of the 
system (6). To this end we look to develop the integral representation for required 
solutions of Laplace integral type 

E(x,y,z) = f / -^exp(in y y + in z z)B(x,n y ,n z ) , (7) 

JC 2,TT J-oo 271 

where E (x,n y ,n z ) is assumed analytic in some domain and the path of integration in 
the n y plane is such that the integrand vanishes rapidly at the ends of the contour C or 
at infinity. Next 



where n opt = Jq /(1 + q ) and 

= cos((p)n y — sm(ip)n z (9) 

rig = sin({p)n y + cos(tp)n z . 

We cannot easily deal with the system (8) as it stands after transformation, but recalling 
the conclusions of ID theory [HE] that the effective conversion is possible if both <C 
0(1) and ~ ±n opt 0(0) we may reduce it to a simpler form. Expanding around 
n opt and around we have found that n° y = sin ((f)n opt and n° z = cos (ip)n opt . Keeping 
in (8) terms being the first order quantity one has 

i opt / Q \ X „ 

where the explicit expressions for Sn^ and Sn^ could be found by varying (9). Making 

the backward transformation with respect to 8n y yields 

_, , / n n \ f°° ddn z 

E(r) = exp [in z + in"y) \ — — exp (ion z z) x 

V 1 J-oo Z7T 

dy'G (y - y') exp (-itan((p)5n z y')E(x,y' , Sn z ), r = (x,y,z) (11) 

G(y-y') = e ^ ( - i(y -f /R2 \ J? = 4L^'|sin M | (12) 
yiixtt 



where the electric field's components E = (E + , eA obey the system of equations 

MT X - ic ° s 4) n °* E+ -t E < =o - (i3) 
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Introducing new notations 

2 1/4 ( y \ d d 

r i /2 1/4 • x, j-^r + 2n opt L by Sn z -> x, y, d± = — ± i— 

Lnql 1 \ cos (v?) / dx dy 

L n q cos ((p) ~ ~ . E + 

a — — , £/+, ^ — > £/+, p — 



^ l + go V 1 + go 

we read (13) in the compact form 

d+Ei. + (x - ay) F = 

d-F -xEt = Q (14) 

In order to make the first step in analysis of the electric field components behavior 
in the conversion region, we consider in the next Section WKB solution to the system 
(14), namely we focus on the quality analysis of the ray trajectory along which the 
energy of the incident beam of the ordinary (or extraordinary) waves is carried over. 
Although in the mode conversion region the WKB approximation breaks down and a 
local full wave equations (14) to be solved, the ray trajectory analysis could be useful. 



3. Ray Hamiltonian dynamics 

The ray representation of the wave field in four dimensional ray phase space r = (x,y), 
n = —id/dr is governed by Hamilton's equations: 

dr dD_ 
ds dn 

where 

D (r, n) = (n 2 x + 8n 2 y - x{x - ay)) /2 = 0, (16) 

is the local dispersion relation and s denotes the orbit parameter. Since the general 
picture of the ray behavior in the two-dimensional subspace xy is of interest, let us 
re-parameterize s — > r, where dr = ds \dD / dn] -1 . Substituting (16) into (15) we obtain 

fr dU(r) 

dr 2 dr ' 1 ' 

where U = x(x — ay)/2. The equations constituting (17) are coupled with a being the 
coupled coefficient. One can introduce the normal co-ordinates r = (x,y) — > (u,v) 
reducing the system of equations (17) to two independent ones. The set of linear 
independent solutions for either of the two is 

w~exp(±z/r), v= (Vl + a 2 + l) 1/2 /V2, 

u~exp(±iwr), tu = (y/l + a 2 - if' 2 /V2. (18) 

The expressions (18), which show the projections of Hamiltonian rays on the (u, v) plane, 
deserves a few comments. First, note the oscillatory behavior along the v direction, 
while the u- motion displays the influence of a retarding force. Second, due to \v\ ^> 
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oscillations along the v direction are not to be clear-cut. Third, we can anticipate the 
full wave solution to the set of equations (14) could be represented as the superposition 
of the eigenmodes intrinsic to the confining potential (~ v 2 ) in the v direction leading to 
oscillatory v motion and the superposition of the linear independent parabolic cylinder 
functions along u direction, which proper for the effective potential ~ — u 2 . The full- 
wave solution to the set of equations (14) is studied in the next Section. 



4. An integral representation of solution 

We would like to obtain the appropriate particular solution of the system (14) matching 
to the WKB solution outside the mode conversion region 

oo 

E| u ^ ±00 ~ eexp (±w 2 / 2 ) H^p(^), (19) 

p=0 

where e is the polarization vector, ip p {v) is a set of eigen functions to be found, 
exp (— iu 2 /2) and exp(w 2 /2) correspond to incident O mode (OX conversion) and 
incident X mode (XO conversion). 

With that end in view we introduce two possible functional substitutions 
providing the identity of the equations constituting (14) 

F = lx±-f J- d+) W ± (20) 

\ y/l — ia / 

/ *4 + s \ r (21) 

\ y/l - la J 
with W + and W~ being required functions that satisfy the equations 

+ + (x 2 - axy ± i/VT^) W± = 0. (22) 

Let us use the transformation 

x = (l — N 2 ^j ^ (cos (ip)u — sin {ip)v ) , 

y = (l - N 2 ) 1/4 (sin (ip)u + cos (V» , (23) 

where cos (tf) = (1 + N 2 )~ 1/2 , sin (ip) = N (1 + N 2 )~ l/2 and N = -(VTTa 2 - I) /a, 
that makes the potential in the brackets in (22) symmetrical 

^ + ^+(* 2 -^ 2 ± (iV + 0)^ = 0. (24) 
If a <, > (see figure 4) then N >, < 0. Seeking the solutions to (24) in the form 



W±{u, v) = J2 W v («) M V (25) 
P =o v J 

where the definition for 4> p (v) is given by 

My/\N\v) = y i/4 exp (-N 2 v 2 /2)H P (y^v>) , (26) 



(2Pp!) 1/2 7T 1 /4 

<P P N\N\v\ <p k (J\N\vj dv = S pk 
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and Hp are Hermitian polynomials, we obtain the equation for and W 

= 0, (27) 

UU~ \ 71 / 

where 

, ir\N\ ( N\ 

Since the properties of 7^ are important for subsequent analysis we list theirs below 
expressing its argument N explicitly 

7p + (N > 0) = 7; (iV < 0) = 7r[AT] (p+ 1) , 

7p (iV > 0) = 7p + (iV < 0) = n\N\p, (28) 
7 P + _ X (JV > 0) = 7p (N > 0) , 7p + +1 (iV < 0) = 7 P " (iV < 0) . 

The solutions and W^ - to (27) matching to the WKB solutions (19) for the 
incident ordinary and extraordinary waves (see Appendix A) are 

W p + (u) = B+D^ (v^exp («r/4)u) (29) 

W~(u) = B p D_ n - /7T (-V2 exp (-i7r/4)«) , (30) 

where B^ are an arbitrary constants, which we choose so that (20) and (21) fit the 
incident WKB wave outside the conversion layer. As it is demonstrated in Appendix A 
W p + and W p ~ correspond to the incident ordinary and extraordinary waves, respectively. 
The parameter 7^ has meaning the length of an evanescence layer (see the equation 
(A. 3)). Its value depends on the mode's number p and combination of the parameter's 
iV sign and the direction of the process (see (28)). That is, first, either of the mode 
tunnels through an evanescence layer with the efficiency inherent itself, second, the 
conversion efficiency being in ID model the invariant under the conversion's direction 
reversal and the poloidal position of the incident beam changing has no longer the same 
property in 2D inhomogeneous plasma. Being mentioned for the first time in jH] for the 
oversimplified plasma model ignoring the poloidal magnetic field and confirmed in [H] 
for the reasonable plasma model accounting for the poloidal magnetic field B y /B 1 
this property depends entirely on the essential two-dimensional character of the waves. 

Inserting (29) and (30) into (20), (21) and (11) after simple algebra we obtain an 
integral representation of solution to (6) 



E+(r) = exp (in° z z + in° y y^j yjl + q x 

" exp (iSn z z) (±<I ± > + <R ± >) (31) 



2tt 

E^(r) = exp (in° z z + m°y) VT+ 
d5n 



la x 



/ — ^ exp (i5n z z) (=f < I ± > + < > 
J-00 2ir ^ ' 

/oo 
dy'G (y - y') exp {-i tan (<p)5n z y') (...) , (32) 
-00 



< ... >= 
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where upper signs at < > correspond to the incident ordinary wave and lower signs 
correspond to the incident extraordinary wave. The explicit expressions for ^ and 
are 



J ± = E \/\N\v)D^ /n (iv^exp (±i7r/4)«), N > 

fc=0 v ' 

OO , s 

I ± = Y,C k <j) k U\N\v)D ±l ^ /7T (±V2exp(±t7r/4)u), N<0 (33) 

fc=0 



and 



i? + = 2|iV|exp ) x 



4, 

oo 

E Ckk<i>k-i(yf\N\v)Dvy-,«-i (v^exp (ztt/4)u) , iV > 

fc=0 

i?+ = -y^exp x 

oo 

E ^ k+ i(#)\/ ir -i (v^exp (i7r/4)u), TV < 0, 

fc=0 

iT = 2|iV|exp x 

oo 

E Ckkfa-^y/lNWD^^ (-v^exp (-i7r/4)«), iV < 



fc=0 

i?- = -^jiV| exp x 

E Ckfa+^y/lNWD^,^ (-v^exp (-i7r/4)«), iV > 0. (34) 



fc=0 

While outside the mode conversion region at u — > ±oo < J + > and < I > match 
to the geometrical optics solutions that correspond to the WKB incident and conversed 
waves (compare the first terms in the r.h.s. of (A.l) and (A. 2), (A. 6) and (A. 7)), the 
WKB asymptotics of < R + > and < R~ > correspond to the WKB wave reflected from 
the evanescence layer (see (A. 5)). 

It is possible to obtain asymptotic expansion of the solution for (p — > 0. For the 

incident ordinary wave at > we have 

, . . /■oo dbn J 30 , 

E+(r) = J I + g exp (in° z z) / — — E ex P (iSn z z) C k x 

V ' J-oo 2tc k=Q 
(M\fW\v)D^- /n (V2exp {m/A)u) + 



2\N\ exp ^JA;0 fc _ 1 (^iv]t;) J D l7 - A _ 1 (v / 2exp («r/4)u)) (35) 
Ee(r) = y/l + ia exp (in° z z) / V" exp (i5n z z) C k x 



k=0 



(-M^\N\v)D, n - /7T (V2exp («r/4)ii) + 



2|7V| exp (i^k^JlN^D.-.^iV^exp (in/4)u)). (36) 
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Last expressions recover the expressions for the electric field components derived in 
for plasma with B y /B <C 1. 

5. The asymptotic expressions for the electric field's components 

For usual spherical tokamak's configuration O mode cut-off layer and the X mode cut-off 
layer cross at the small angle, i.e. a <C 1. Keeping only terms in (23) being zero and 
first order quantity with respect to the parameter a we read it as 

x — u — av/2, 

y^v + au/2, (37) 
Then we simplify 

oo 

< i± >~ Yl ••• < <t>k{y/d) >, 

fe=0 

oo 

< r± >~ Yl ••• < (y/ d ) >> 

fe=0 

where d = (2 1//4 (1 + g ) 1//2 L^ 2 cos(t^) 1 / 2 ) / {iOnzX} 1 / 2 ^^ and y is scaled in the units 
A = 2i\cjuj. Comparing the asymptotic representation of D±i 7p / n for the argument 
x — > oo and asymptotic representation of D± ilp / n for the argument x — > — oo (see 
Appendix A) we obtain an explicit expressions for the conversion coefficient, Tox,xo, ( m 
energy) and reflection coefficient, Rox,xo, where the subscribes OX and XO correspond 
to the direction of the process 

1 °° roo fj^n , x 

Tox=pE ^IC/exp (-2 7 J , Rox = l- Tox, (38) 

p=0 J-oo ZTT 



1 x roc rjm . . 

T xo = jsY, ^\C p \ 2 e W (-2^),R xo = l-T xo . (39) 



and 

The coefficient C n is given by the expression 
C p = (2tt) 2 x 

/ A(y, z) exp (-i27f(n° z z + 5n z z + n Q yj) < <f)* p {y/d) > dydz, (40) 

J —oo * " ' 

where * means the complex conjugation, A(y, z) is given field distribution in the incident 
beam in the WKB region, x and z are scaled in A and P is a normalizing coefficient 
proportional to the beam intensity. The expressions (38) and (39) derived similar to 
ones obtained in (Hj El excepting the formula for C p given by (40) and deserve one 
comment. Noting Tox ^ Txo at the fixed a, we could read the reciprocal relation for 
the conversion coefficients in 2D inhomogeneous plasma as 

Tox (B) = T xo (-B) , R ox (B) = R XO (-B) , (41) 
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where the argument B is presented explicitly. For the first time the property (41) 
was mentioned in |BJ. In jU] it was proved for plasma with poloidal magnetic field 
B p /B < 1. 

For illustration, in figure 5 the conversion coefficient, Tqx, of an incident beam 
with transversal distribution 

A(y, z) = (vrp 2 ) ~ V4 exp (-y 2 /{2p 2 ) + in opt sin (<p)y + in opt cos (ip)z) , 

calculated with the use of ID model [4 (dashed curve) and the formula (38) (solid 
curves) for a < and different <p versus the dimensionless radius, p y /\, is shown. 
Calculations are executed under conditions of EBW heating experiment in MAST 
L n /\ = 3, L by /\ = 70, q = 0.7, u/(2tt) = 15GHz. Rather strong variation of 
dependence presented in figure 5 with variation of (p is observed. 

6. Summary and conclusion 

The mode conversion of O mode to X mode have been examined in a spherical tokamak 
geometry for plasma with a cold plasma dielectic tensor. We have assumed the 
wavelengths are much shorter than the equilibrium plasma gradient length, L n . Thus, 
we expect WKB approximation to apply except possibly near O-mode cut-off surface 
or mode conversion layer. 

Expanding wave equations in the region near the intersection of O and X mode cut- 
off surfaces, one finds reduced wave equations appropriate for this region which depends 
on two coordinates, x and the poloidal angle y. We seek the required solution to the 
set of equations so that it matches to the WKB solution outside the mode conversion 
region. After functional substitution we have found the integral representations of the 
required solution. 

To the case typical of a spherical tokamak's configuration for which the O mode cut- 
off layer and the X mode cut-off layer cross at the small angle we have been simplified 
the obtained integral representations of solution and found the wave fields and the 
conversion coefficients explicitly. 

We have shown that at the fixed poloidal position, i.e. fixed a, the conversion 
coefficients obey the reciprocial relation (41) as it was shown earlier in the frame of the 
model neglecting the poloidal magnetic field [H] and of the model in which the poloidal 
magnetic field is assumed to be small, B y j B <C 1. 

The importance of 2D effects has been demonstrated by the example under the 
MAST conditions. 
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0,4 




Figure 1. ip = arctan Bg/B z behavior versus the poloidal angle for a typical MAST 
tokamak discharge 




B 



Figure 2. Beam with radius p of ordinary polarized electromagnetic waves incident 
on the O-mode cut-off surface. The co-ordinate ro indicates the position of the beam 
center on the O-mode cut-off surface 




(a) (b) 



Figure 3. Illustration of the co-ordinate systems (x,y,z) and (x, C, £) used, (a) Two- 
dimensional subspace (xoy); (b) Two-dimensional subspace (yoz) with ip being an 
angle between the magnetic field and the toroidal direction 
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(a) (b) 



Figure 4. Co-ordinate systems (x,y) and (u,v) (a) a > (N < 0); (b) (a) a < 
(N > 0) 




p (0/ (27Ec) 

Figure 5. The conversion coefficient, Tox, versus the dimensionless radius, p y /X, 
of an incident beam calculated with the use of ID model [3] (dashed curve) and the 
formula (38) (solid curves) for a < and different ip under conditions of EBW heating 
experiment in MAST L n /X = 3, L by /X = 70, q = 0.7, w/(2tt) = IhGHz. 

Appendix A. Asymptotic expressions for the parabolic cylinder functions 

We pay our attention to the functions of the parabolic cylinder with different arguments 
and indices D^/^ (in /4)uj and D^/^ y^exp (— m/A)uj. Let us begin from 

asymptotic expressions for the first one at two limits u — > ±oo which describe the 
corresponding WKB solutions beyond the conversion layer. When the argument u tends 
to +oo [Hj 

A 7 /tt (V2 exp (i7r/4)«) | u ^oo « (v^u)* 77 " exp (- 7 /4 - iu 2 /2) , (A.l) 
while at u — > — oo 

D il/n (v^exp (i7r/4)u) U^-oc ~ (v^wj^exp (37/4 - iu 2 /2) - 
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\J~2~k exp(z37r/4) 



r(- 



^AO (V2u) 1+n/7T 



exp 



(j/4 + iu 2 /2). (A.2) 



The second term in (A.2) is much smaller (with factor l/u) compared to the first one. 
Therefore the function D^/^ (^y/2exp (m/A)uj in the WKB sense describes the wave 
approaching an evanescence layer from — oo and conversed wave propagating from an 
evanescence layer to oo. Comparing the coefficients at exp (— iu 2 /2) in (A.l) and (A.2) 
one obtains the conversion coefficient (in energy) for the wave described by the parabolic 
cylinder equation 



T = exp(-2 7 ). (A.3) 

Farther we obtain the derivative of the parabolic cylinder function 
d 

Q^ D h/n (v^exp (in/ 4) uj = 

— ^exp (z7r/4)Dj 7 / 7r _ 1 [v^exp {iix/A)u \ — iuD^/^ fy^exp (iir/4)u) (A. 4) 



7T 



The first term in the r.h.s. of (A. 4), which has the asymptotics 
A7/71--1 (v^exp (i7r/4)wj | « 0(l/|u|) 
A 7 /7r-i (v^exp {m/A)u \ | u _>-oo ~ 



^7T 



^— ^ exp (7/4 + iu 2 /2) + 0(1/|«|), (A.5) 



r(i-<?) (V2 U y 

describes the reflected wave. 

Comparing the asymptotic expansions for the parabolic cylinder function 
£>_ i7/7r (-v^exp {-m/4)u) [El 

D_ ij/n [-V2 exp (-Z7r/4)u) l^oo w (vW) 7/ exp (37/4) exp (m 2 /2) - 
\f2~K exp (— i37r/4) 



r(i7/7r) 



^exp( 7 /4)exp(-m 2 /2). (A.6) 



and 



D_ il/n (-v^exp (-in/A)uj |„_-oc ~ 
v^m) exp (-7/4) exp (iu 2 /2) , (A.7) 



we note that it describes the wave which incident on an evanescent layer from u — > 00 
(see the first term in (A.5)) and propagates from an evanescent layer to the plasma 
boundary {u — > —00) (see (A.6)). 
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